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SYMBOLS 


p 

beam waist diameter (to 1/e ) 

product of fringe number and time, C = N e ff* t g ^ 

skin friction coefficient, = x/q e 

spot location parameter, C x = X 1 /X 2 

distance from focal plane to surface 

error 

lens focal length 
oil film thickness 
beam intensity 
index of refraction 

fringe number; number of optical path lengths in the oil film 
photodiode output, proportional to light intensity 
dynamic pressure 
fringe visibility 

beam separation due to beamsplitting window 
time 

velocity components in Cartesian system 
distance from beam waist to lens 
Cartesian coordinates 
oil-film wedge angle 
thickness 

optical path difference due to beamsplitter 


wavelength 



viscosity 


li 

8 beam incidence angle (measured with respect to normal local surfaces) 

p density 

x skin friction or shear at air-oil interface 

Subscripts 

A,B,C,D,E pertains to light paths depicted in figure 1(b) 
a property of the air 

e pertaining to the freestream of the air flow 

eff "effective" value 

l characteristic of the laser light source 

in input beam properties 

0 property of the oil 

out output beam properties 

w property of. beamsplitter window 

<f> value at oil-film leading edge or initial value 

1 "upstream"-most measurement position 

2 "downstream"-most measurement position 
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IMPROVED SKIN FRICTION INTERFEROMETER 


R. V. Westphal, W. D. Bachalo,* and M. H. Houser* 
Ames Research Center 


SUMMARY 


An improved system for measuring aerodynamic skin friction which uses a dual- 
laser-beam oil-film interferometer has been developed. Improvements in the optical 
hardware provided equal signal characteristics for each beam and reduced the cost 
and complexity of the system by replacing polarization rotation by a mirrored prism 
for separation of the two signals. An automated, objective, data-reduction proce- 
dure was implemented to eliminate tedious manual manipulation of the interferometry 
data records. The present system was intended for use in two-dimensional, incom- 
pressible flows over a smooth, level surface without pressure gradient, but the 
improvements discussed are not limited to this application. 


INTRODUCTION 


A thin film of oil on a smooth surface will flow under the action of the aero- 
dynamic shear stress created by the flow of air over the surface. The equations 
governing the behavior of an oil film in laminar flow are well known (refs. 1-3); 
figure 1(a) depicts a thin, two-dimensional oil film flowing under the action of 
aerodynamic skin friction. The skin friction is assumed to be steady in time and 
constant in the stream and span directions. Skin friction gradients and pressure 
gradients (including gravity) as well as surface curvature effects are not included; 
these effects can be included if they are present (refs. 1 and 2). A simple rela- 
tion then exists between the time-dependent oil film thickness h at a location 
X and the skin friction: 


t = [(X - X (() )y o ]/[h(t,X)t eff ] (1) 

As this equation gives the time-dependent variation of oil film thickness at a known 
location and the oil viscosity is known, the skin friction may be calculated from 
equation ( 1 ) . 

A number of questions remained from previous investigations regarding the 
response of a flowing oil film to skin friction produced by a turbulent boundary 
layer. (1) How do the fluctuations in skin friction which are encountered in real 
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turbulent flows influence the oil film? (2) Why are "waves” observed in oil 
films? Both of these basic questions were examined in detail in a companion numeri- 
cal study (ref. 3) in which it was concluded that the oil film will produce a cor- 
rect measurement of average skin friction if the average rate of thinning is mea- 
sured, and that the "waves" are formed from the interaction of the oil film with the 
turbulent boundary layer. The study also showed that near the oil-film leading 
edge, there must always be a wave-free region where no interaction takes place and 
the linear lubrication theory (eq. (1)) is valid. This has also been observed 
experimentally, although in supersonic flows Monson has found that this region 
becomes vanishingly small. 

To put the method into use for quantitative measurement of skin friction, a 
technique for accurate film-thickness measurement is required. Wave-free flow at 
very low Reynolds numbers and the need to confine the film to an area of a few 
centimeters require oil-film thicknesses of less than 50 pm. Under these condi- 
tions, optical measurement using the principle of interferometry is the method of 
choice. 

Although the principle of interferometric film thickness measurement has been 
known for some time, e.g., Jenkins and White's text (ref. 4), Tanner and Blows 
(ref. 2) at Brighton Polytechnic were apparently the first to examine quantitative 
interferometric measurement of oil film thickness as a means for direct determina- 
tion of skin friction. Murphy and Smith (ref. 5) attempted to apply interferometric 
film-thickness measurements to obtain evaporation rate data which were to be related 
to skin friction through a simple mass-transfer analogy. Later, Wazzan et al. 

(ref. 6) attempted to apply the same technique to a supersonic flow. They actually 
worked out the relationships of lubrication theory needed to relate the rate of 
thinning of a liquid film to skin friction, but appear to have overlooked the use- 
fulness of the idea. 

The principle of interferometry with a laser beam can be used to measure the 
time-dependent oil film thickness as shown in figure 1(b). Some of the incident 
beam is specularly reflected at the air-oil interface (ray A); the remainder passes 
through the oil film and is either transmitted through (ray B) or is reflected by 
the solid surface (ray C). The reflected light (ray C) passes through the oil and 
is partially reflected again at the oil-air interface (ray D — an internal reflec- 
tion). Multiple internal reflections can be ignored for small reflectance 
(ref. 4). The transmitted portion of the primary beam (ray E) will interfere with 
the first (specular) reflection from the air-oil interface, thereby providing a 
light signal. As the phase relation between A and E changes, the intensity of the 
signal will vary approximately sinusoidally because of a variation in the optical 
path length. This variation is caused by a thinning of the oil film. One 
wavelength in the oil is X Q (one "fringe"), corresponding to a thickness variation 
given by references 2 and 4. 


X = X / [ 2n cos 9 ] 
o l o o 


( 2 ) 
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where 


e = sin _1 [(sin 9„)/n ] (3) 

o i o 

Let N denote the film thickness in number of optical fringes, then X Q N(t,X) is 
equal to the oil film thickness h(t,X). Equations (1), (2), and (3) give the basic 
equations needed to determine skin friction once the time-dependent fringe history 
(thickness) of an oil film has been measured at a particular location X. Note that 
two measurements a known distance apart will be required if the precise location of 
the oil-film leading edge is not known. 

The signal-to-noise ratio is strongly influenced by fringe visibility r, 
defined as the ratio of the intensities of the two interfering beams 



Obviously, maximum visibility is obtained for equal intensity of the interfering 
beams. The fringe visibility is influenced by the optical properties of both the 
oil and the surface, as is further discussed later. 

Tanner's instrument (refs. 2,7-12) consisted of a fairly small (less than 
50 mW) He-Ne laser, beam steering and focusing optics, and receiving optics 
(fig. 2). A beamsplitter was used to give a second "reference" spot for locating 
the measuring spot with respect to the oil leading edge. Photodetector output was 
monitored on a strip-chart recorder and later timed using an analog fringe-detection 
circuit and a computer (ref. 13). Very long runs were taken so as to thin the oil 
to less than one fringe. Then, the time at which N = 1 occurred could be deter- 
mined, and a reference fringe number was obtained to allow accurate fringe number- 
ing. Alternatively, the origin for time was taken as the start of the run and the 
reference fringe number was chosen to give a nearly constant fringe-time product. 
Measurements were made on coated and uncoated glass as well as on polished metal 
surfaces. Tanner (ref. 13) also evaluated techniques for viscosity measurement 
which can be used when oil viscosity is not accurately known. Tanner and his 
co-workers have done extensive comparisons of the oil-film interferometer technique 
with the classic Preston tube (refs. 9 and 10), and have also used oil film inter- 
ferometry for visualization of surface flow patterns (e.g., Tanner, ref. 12). 

A "surface phenomenon" was investigated by Tanner (ref. 13). The results indi- 
cated that extremely thin films of silicone oil may not accurately follow lubrica- 
tion theory. These results imply (1) that small-surface roughness (0.2 ym or less) 
can give rise to an apparent "stagnant" region near the surface and (2) that the 
oil-film leading edge position is best represented by an "effective" value — probably 
the result of surface tension. Potential difficulties stemming from (1) are circum- 
vented by the use of very smooth polished surfaces, and by terminating the measure- 
ments with a film thickness of approximately 1 ym or so. For the dual-beam system, 
(2) is automatically satisfied. 
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Monson and others at NASA Ames Research Center (refs. 14-17) later extended the 
method by introducing a second data beam, which eliminated the need to locate the 
leading edge of the oil film. The two data beams were perpendicularly polarized and 
were separated using polarizing filters. Otherwise, the optical system that was 
used (fig. 3(a)) was similar to Tanner's. Monson 's data-acquisition procedure was 
to record the output of each photodetector on separate channels of a strip-chart 
recorder for a period of time sufficient to produce about 20 "good" fringes from the 
downstream photodetector. The times elapsed for 10 and 20 fringes on the downstream 
signal were used to calculate the effective time and fringe numbers, thus removing 
the requirement of using lengthy data records to provide a reference fringe number 
or of measuring or estimating the time origin. The time elapsed for 10 fringes at 
the upstream location gave the fringe reference for the upstream beam. The results 
for the two beams were then combined to give the skin friction. The measurements 
were made on polished metal test surfaces (aluminum and stainless steel). "Bad" 
records, defined as those which. contained uneven visibility or obvious noisiness, 
were identified qualitatively by direct inspection of the strip-chart recordings. 
Figure 3(b) shows a typical signal from Monson 's instrument along with the data- 
reduction equations that were employed for two-dimensional flows with no pressure or 
shear gradients. They have also used the instrument in a wise variety of flow 
conditions such as three-dimensional boundary layers, separated flows, and complex 
supersonic flows. 

The present application to skin friction measurements in low-speed, two- 
dimensional turbulent boundary layers with zero pressure gradient required a mea- 
surement uncertainty of ± 5 %. Discussions with Monson (private communication) and 
the authors' initial tests of a prototype single-beam system revealed a number of 
problems with the existing instrument. Usual data-reduction procedures called for 
tedious and somewhat subjective manual manipulation of fringe records. "Bad" runs 
were often obtained and later discarded, and the fringe visibility index was quite 
low (r < 0.2) because metal test surfaces were used. Monson 's system gave a differ- 
ent signal quality for each beam (fig. 3(a)), which tended to complicate processing 
of the fringe records. Monson also indicated that oil heating could be produced by 
the use of highly focused beams if too much laser power was employed. 

Two major objectives evolved for the current work. The first was to design and 
build a new two-beam interferometer that would give matched signal quality on both 
beams. Maximum fringe visibility was sought to give the best signal-to-noise ratio 
and to allow the' use of the least amount of laser power (thereby minimizing oil 
heating). The second objective was to automate the fringe processing procedure to 
remove any subjectivity in the identification of bad records and to speed the data 
processing so that skin friction could be measured on-line during an experiment. 
Initially, an uncertainty analysis was carried out to establish requirements for 
tolerances and design parameters. This analysis is described in the next section, 
followed by sections which outline the optical hardware and data-reduction proce- 
dures. Problems encountered with the system, as well as sample results, are given 
in the Discussion section. The Concluding Remarks include the authors' plans to 
implement further refinements to the interferometer. 
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UNCERTAINTY ANALYSIS 


An uncertainty analysis was undertaken to aid in the design and to point out 
critical measurands. Combining (1) and (2) and then introducing the (constant) 
product of fringe number and time C gives 

T. = v (X - X.)/[x C(X)] (5) 

O 9 O 

Applying (5) for each beam to eliminate X^ gives the dual-beam data reduction 
equation 


t = u o (x 2 - x 1 )/{x 0 [c(x 2 ) - C(X 1 )]} (6) 

The uncertainty in t can be evaluated by a simple root-sum-square combination of 
the uncertainties from the individual terms under the- assumption that these terms 
are mutually independent (ref. 18). For the present application, the uncertainty is 
expressed as a percent of the indicated value and all values are estimated at 20:1 
odds. The analysis is summarized in table 1, in which the partial derivatives and 
uncertainties in each term are individually evaluated. Four separate terms are 
identified and discussed below. 

Oil viscosity appears as a linear term in equation (6) and therefore contrib- 
utes strongly to the uncertainty. Tanner (ref. 13) and Monson (ref. 14) both 
directly measured the viscosity of the oils they used; Monson' s measurements were 
used in this analysis. Fifty-centistoke (cs) (nominal) oil such as Monson used 
(ref. 14) was obtained, the viscosity of which was found to conform to within 2 % of 
the manufacturer's quoted values (see fig. 4 of ref. 14). The manufacturer's rela- 
tion for viscosity was adopted 

“o = “nom «P<-0.01M<T o - T nom )) (7) 

Oil temperature was assumed to be equal to the flow temperature, which was continu- 
ously monitored during each experiment. When the temperature drifts or the flow is 
not isothermal,, special treatment is required as the oil viscosity changes rapidly 
with temperature. Generally it is necessary to measure the oil temperature more 
directly (e.g.„ with a temperature probe embedded in the surface below the oil). 

Beam spacing (X 2 - X^) also appears as a linear term in equation (6), and was 
measured directly to an estimated uncertainty of 0.05 mm, or about 2 % for the small- 
est spacing used (2.77 mm). The measurement was made using a diffuse surface placed 
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at the spot focal plane. Repeated trials were averaged to yield the quoted 
uncertainty. 

Alternatively, if oil viscosity was measured directly by one of the techniques 
employed by Tanner (ref. 13), the instrument could simply be calibrated in terms of 
the product of beam spacing and viscosity. This procedure may somewhat reduce the 
net uncertainty. 

The uncertainty in the optical path length in the oil depends primarily upon 
accurate knowledge of the oil index of refraction, because the laser wavelength is 
known precisely and the incidence-angle term cos 0 Q is very nearly unity for near- 
normal incidence. The oil index of refraction, which was nominally 1.4, was mea- 
sured as 1.405 using an Abbe refractometer . The uncertainty in this measurement is 
estimated to be less than 0.5# — leading to an uncertainty in X Q that is negligible 
compared to the other terms. 

The fringe-time products C(X ^ ) and C(X 2 ) are the primary measurands of the 
system. The difference of these two quantities has been considered a single term 
for convenience. It is apparent from table 1 that the net-uncertainty contribution 
from this term depends upon the selection of measurement location if the uncertainty 
in the determination of C is taken as a fixed percentage of the value. The 
results showed that as long as more than 20 fringes were used, the value of C 
could be determined with an estimated uncertainty of about 0.5#. This estimate was 
obtained from experiments with two-dimensional flat-plate boundary layers with 
minimal gradients of pressure and skin friction at shear levels of 0.5 to 2 N/m 2 . 

For a 3# to 4# uncertainty contribution from this term, this estimate implies that 
the value of (X^ - X^) should be less than about four times the beam spacing. 

The minimum value of (X 1 - X^) was determined by the requirement that more than 
20 fringes be identified. The limited spatial resolution caused by the finite spot 
size of the beams caused reduced fringe visibility if the oil-film wedge angle was 
greater than approximately 0.1°. Twenty fringes correspond to a thickness of 
approximately 5 um. This thickness implies that the minimum distance to the leading 
edge should be approximately 3 mm, which was approximately equal to the smallest 
beam spacing that was used. 

These uncertainties indicate that the minimum uncertainty likely to be attained 
with the oil-film-interferometry technique is about 3#, which is primarily caused by 
the uncertainty in the determination of the oil viscosity. The analysis shows that 
the main contribution to the uncertainty for the current instrument is the uncer- 
tainty in the determination of the difference in the finge-time product at the two 
measurement locations. This uncertainty may be limited to 3# to 4# if 
(X 2 - X^) < (X^ - X^) < 4(X 2 - X^) is maintained. The overall uncertainty in 
measurement of skin friction, subject to the assumptions outlined above and valid 
for the system described below, is estimated at approximately 5#. 
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OPTICAL HARDWARE 


A photograph and schematic of the optical system are shown in figure 4. As in 
previous configurations, the current system transmits two beams focused onto the 
test surface and uses PiN photodiodes as the detectors. Both transmitter and 
receiver components share a common optical-rail assembly. The test beams pass 
through a slot in the rail with the test surface located below the optical assem- 
bly. Light scattered from the surface is collected through the same slot so that 
the test beam incidence must be near-normal to the surface (actual incidence angles 
were about 3°). The key optical-design change from previous work is the physical 
separation of light that is reflected from the two surface spots through the use of 
imaging lenses and a reflective edge. This change allowed for simplification of the 
optical system and maintained good separation (minimal crosstalk) between the two 
channels. 

A low-power (0.8 mW), polarized, helium-neon laser provided the coherent light 
source. The beam is directed by mirrors through a beam-expanding and -focusing 
telescope, composed of a 20-power microscope objective and a 150-mm-output lens. By 
adjusting the separation between the objective and output lenses, a focused spot was 
produced on the test surface at a considerable distance from the output lens. The 
equations which describe the transformation of the (approximately) Gaussian laser 
beam are 


fa. 

in 


out 2 21/2 

« ira. 

f)2 t _i" 


( 8 ) 


(x. 

in 


out 


= f + 


(x in - f)f* 


2 2 

, p .2 wa in 

x in ' f) + X“ 


(9) 


Applications of these equations with values for the actual optical configuration 
yielded a focused beam diameter of approximately 42 pm. 

After leaving the beam expander, the beam is directed onto a plane-parallel to 
10 arc-sec) window of known thickness. At an incidence angle of 45°, both the front 
and back surfaces of the window reflect about 4$ of the incident beam. This 
arrangement splits the expanded beam into two test beams of nearly identical inten- 
sity, with a separation given by the following equation: 
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( 10 ) 


26 sin 9 . 
w % 


tan cos 


-i 003 \ 


A fixture was fabricated that allowed windows of different thicknesses to be inter- 
changed, which produced spot separations in the range of 2 to 10 mm. 

The beamsplitting method just described introduced a difference in optical path 
length between the two test beams. This difference, which changes with the window 
thickness, is given by 


A 


2n 6 , 
w w 


sin cos 


-i 003 \ 


w 


s 

tan e„ 

a 


( 11 ) 


The consequence of the path difference is a change in focus location, which prevents 
focusing both beams simultaneously on the test surface. It is not critical that 
focus occurs at the test surface, but it was considered desirable to produce two 
spots of equal size. In practice, the beam expander was adjusted to equalize the 
resulting spot sizes at the test surface. This "out-of-focus" condition produced 
one beam focused above the test surface and the other beam focused below the test 
surface. The change in beam diameter with distance from the focal plane is given by 


a(d) 


a . 
in 


Xd 


2 1/2 


ira . 


m 


( 12 ) 


Most of the testing was done with a BK-7 window of a nominal thickness of 3.2 mm, 
which introduced a path difference of about 8.5 mm. Proper adjustment of the beam 
expander would therefore produce two spots on the surface, each located 4.25 mm from 
the focal plane. The results of the application of equation (12) indicate that this 
condition will produce surface spot sizes of 92 ym — about double the focused spot 
size. 


After the beamsplitting was performed and the component losses were taken into 
account, each of the two test beams contained less than 3 % of the original laser 
power — 24 yW per beam for the present system. This resulted in a surface-power 
density of less than 0.004 W/mm 2 per spot, which minimized oil-film heating but 
still provided sufficient signal amplitude for the detectors. 


The receiving system consisted of a 450-mm focal length lens which collimated 
both specular and diffuse reflected light. A second lens (f = 250 mm) focused the 
light from the two spots onto opposite sides of the edge of a reflective-coated, 
right-angle prism. Sufficient image separation and spot focus allowed redirection 
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of the light from the two surface spots to individual detectors. An adjustable iris 
and laser line filter preceded each detector to minimize background light detec- 
tion. The photodetectors were identical photodiodes of the PiN type with active 
areas of 1 cm^ and sensitivity of 0.4 yA/yW. The photodiodes were operated in the 
photovoltaic mode by identical, battery-powered amplifiers and had a saturation 
level of 50 mW/cm^. An amplifier gain of 10 V/yA was used, which yielded signal 
levels of a few volts that were convenient for direct input to the computer analog- 
to-digital (A/D) converter. 


SIGNAL PROCESSING AND DATA REDUCTION 


The computerized data-acquisition system (shown schematically in fig. 5) con- 
sisted of a PDP 11/44 minicomputer with 20 Mbytes of hard disk storage and a clock- 
driven, fully differential eight-channel A/D converter. The analog outputs from 
each photodiode amplifier were simultaneously read directly by the A/D converter for 
a period sufficient to record 4000 samples per channel. The sample rate was 
selected (through variation of the clock frequency) to record a sufficiently long 
fringe record so that a thinning of the oil to about 1 ym could be achieved at the 
measurement location farthest downstream. These data were written directly to hard 
disk, then processed on-line with a separate data-reduction program that included 
graphics capability. 

Figure 6(a) shows a typical raw digitized data record from the two photodetec- 
tors. The data-reduction procedure described in this section details the method 
used to process such records to obtain the value of C (the constant product of 
fringe number and time discussed in the section on uncertainty analysis). Note 
particularly that each channel is treated independently, i.e., the product of the 
fringe number and the flowing time for the upstream and downstream channels are 
computed independently by using the same reduction parameters. An unexpected side 
benefit of this fitting procedure is that the fitting parameters give an objective 
measure of the data quality and therefore a criterion for the rejection of bad 
records, as explained further in the next section. Then the results from the two 
beams are combined, along with the optical system parameters, oil properties, and 
tunnel conditions to compute the skin friction. The Appendix contains listings of 
the computer routines (referenced in this section) which were used to perform the 
data reduction,, Other aspects of the data acquisition, graphics, and database 
management software used with our particular hardware are described in a separate 
report (ref. 19). 

In the data-reduction procedure, the fringe locations are identified and then 
the locations (in time) of successive fringes are fitted to the lubrication theory 
with a best-fit criterion. First, the digitized photodetector output is smoothed 
using either an exponential or multipoint smoothing algorithm (both algorithms were 
included in subroutine SCSM00). It was soon concluded that smoothing was not 
necessary because the fringe identification algorithm provided some smoothing and 
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the photodetector output was virtually noise-free. Whenever smoothing was used, 
exactly the same algorithm and parameters were used for each channel. 

A simple fringe-identification algorithm was developed and implemented in 
subroutine SCZERO. A fringe location is most accurately defined as a location where 
the signal amplitude equals its mean level (a "zero crossing," rather from a signal 
peak). This definition is the same as that used by Tanner, but different than that 
of Monson, who used signal peaks. It was decided to use only positively sloped zero 
crossings so that successive zero crossings corresponded to reductions in film 
thickness of precisely one fringe. The zero finder algorithm worked by first deter- 
mining the mean signal level and standard deviation (fig. 6(b)). Next, regions 
where the data values fell within a specified fraction of 1 standard deviation (SD) 
of the mean line were bracketed. Finally, the data in each identified region were 
fitted to a straight line whose intercept with the mean line was taken as the fringe 
location. Note that the algorithm does not account for possible drift in the mean 
line. 


A minimum-error fitting procedure was used to select the best value of the 
fringe-time constant C for each beam (subroutine SCCFIT). Each successive zero 
crossing corresponds to a reduction in oil-film thickness of one fringe, so the film 
thickness in fringes at the kth zero crossing is simply 


W k) * "♦ - k 

Similarly, the effective time of flowing can be defined as 


( 13 ) 



+ t. 


( 14 ) 


where N^ is the reference fringe number and t^ is the time constant. The values 
of N, and t^ were chosen to minimize the relative error (e) in the resultant value 
of C(X), that is, t^N^ so that 


is minimum, with 


e 



( 15 ) 


C k = N eff' fc eff 


and 


max 


C = 


k , , k 
max k=1 


CT = 


1 


max 


max k=1 


10 



Equation (15), with the minimization constraint, is the fitting equation mentioned 
below. 

Both data channels were analyzed by using equation (15), with the minimization 
constraint, which gave reference values and t^ for both data beams. This 
analysis contrasts with Monson's procedure, which computed the effective oil-flow 
time reference from the downstream beam only, and that of Tanner, who generally used 
an actual time of flow. The new method gives the effective time from the least- 
square fit for each channel. The two values of t^ were then compared; substantial 
differences in these constants are taken to indicate bad data, as explained further 
in the next section. Figure 6(c) shows the results of fitting the fringe data of 
figure 6(b) in this manner; it can be seen that 1/ N e ff is precisely a linear 
function of t (; f.f. 

Note that the fringes whose amplitude is less than the specified fraction of an 
SD that is used to bracket the zero-crossing region cannot be identified. This 
means that for records which contain highly variable fringe visibility, fringes will 
be missed. These records were found to yield a poor fit to lubrication theory, 
i.e., a large Value of the error e defined in equation (15). It was found that 
most records yielded e < 0.003, and that e > 0.004 signaled that fringes were 
missed, or that dirt was present on the oil. 

This fitting procedure is applicable only to the case of two-dimensional flow 
with no pressure or skin-friction gradients in the streamwise direction and no body 
forces (such as grvity). If such effects are present, the fitting equation must be 
modified to account for them. However, the data-acquisition and fringe- 
identification procedures just described should still be applicable. A suitable 
minimum-error fitting criterion is substituted for that defined by equation (15). 

For example, appropriate equations are given by Tanner (refs. 2 and 7) for several 
combinations of effects; Monson (ref. 17) has investigated applications to three- 
dimensional flows in detail. 

After C(X^) and C^) are found through the fitting procedure just described, 
subroutine SCLOFI is used to implement equation (6) with the additional input values 
of the free-stream dynamic pressure of the tunnel, the oil temperature, the beam 
spacing, and the incidence angle. The Dow-Corning viscosity relation (eq. (7)) is 
implemented in this routine. 


DISCUSSION AND SAMPLE RESULTS 


As with any measurement procedure, experience in the laboratory with this 
procedure disclosed several unanticipated characteristics. These are discussed with 
the goal of anticipating similar problems in future applications of the method. 

None of these characteristics represented major impediments, but they did affect the 
procedures and equipment finally settled upon. These problems will first be dis- 
cussed; example measurements will then be presented. 
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Method of Oil Application 


Oil is first applied to the surface. Then the tunnel is turned on and allowed 
to achieve the test speed. During this start-up period, the film thickness rapidly 
attains the linear distribution of lubrication theory. For two-dimensional flow, 
the initial, streamwise, oil profile was found through direct calculation to have no 
effect on the results (ref. 3). It has no effect because the linear thickness 
distribution near the leading edge is very rapidly established and the early-run 
data (while the film is relatively thick) are not used. The procedure presented in 
this paper in which the effective oil-flow time and reference fringe number are 
computed from only a segment of the fringe record further contributes to the insen- 
sitivity of the results of the fit to initial conditions at tunnel start-up. 

Theoretically, in the absence of gravity and surface tensile forces, the ini- 
tial spanwise extent and uniformity of the film should also be unimportant. Monson 
(ref. 16) observes that measurements using drops of oil will suffice for measure- 
ments in two-dimensional flows. It was assumed that the presence of both finite 
surface tension and gravity effects would imply the use of oil films that were at 
least several times the beam spacing in spanwise extent. Monson (refs. 14-17) used 
a plastic straight edge which was coated with oil and then pressed onto the test 
surface; this method was found to be tedious in application. To aid in the applica- 
tion of a uniform, initial oil line, a fixture was used (fig. 7) which consisted of 
a syringe and a metal wedge with a small slit. Using this fixture, repeatable and 
uniform oil lines about 4 cm long which contained approximately 0.05 cc of oil could 
be applied. 


Dirt on the Film 

The presence of dirt or dust specks on the oil surface is the single most 
important consideration in the application of the oil film method. Small particles 
can deflect the laser beam or can cause local distortion of the oil flow and conse- 
quent behavior which deviates from lubrication theory. In the experimental design, 
it is desirable to minimize the occurrence of dirt on the oil film surface and to 
identify or note the presence of dirt so that contaminated data records can be 
discarded. Dirt on the film was minimized by using filtered air in the wind tunnel 
and by carefully cleaning the test surface. Use of the oil applicator described 
earlier also helped maintain oil purity by preventing any contact of the dust-laden 
air with the (clean) oil. 

Dirt recognition proved to be a more serious problem. Initially, the oil film 
was inspected visually after each run for specks of dirt. However, about half the 
runs revealed some amount of dirt, although it often was far from the measurement 
location. The dirt-contaminated records always displayed erratic fringe patterns. 
Every attempt to fit these records to lubrication theory resulted in large errors 
and/or unequal values of the time constant t^ from the two beams. This result led 
to the adoption of two criteria for the rejection of "dirty" signal records: 


12 



1. If e > 0.004 for either beam; or 

2 - if is.t - s,2> > °- 02 • ')• 

These relations represent a completely objective criterion for rejecting bad 
data records in a systematic way. Use of the above criteria resulted in a rejection 
rate of approximately 1555 to 20$, which is a significant improvement over the 50$ 
rejection rate obtained from visual inspection of the oil film. If the rejection 
were based on visual inspection of the fringe records for anomalous behavior, up to 
twice as many runs would be rejected if one were discriminating in the inspection. 


Test Surface Finish 

Tanner (ref. 12) has reported that the best fringe visibility occurs if par- 
tially reflective, metal-coated- glass surfaces are used because the oil-free surface 
reflects only 4$ (approximately) of the incident beam. To achieve maximum visibil- 
ity (r = 1), the reflectance at the oil-wall interface should also be about 4$. The 
reflection from a polished metal surface is composed of specular and diffuse compo- 
nents. Surfaces that are highly polished will reflect specularly a higher propor- 
tion of the incident beam. Tanner (ref. 12) reports that excessively polished 
surfaces will give too much specular reflection, thereby causing multiple internal 
reflections which degrade signal quality (see also ref. 4). Only the specular 
reflection from the polished surface will interfere with the (specular) reflection 
at the air-oil interface. It would therefore seem that the polish need merely be 
adjusted to give the proper proportion of specular reflection. 

Polished aluminum inserts were initially used for the measurement surface 
beneath the oil film. However, it proved extremely difficult to obtain a consis- 
tently uniform finish over the entire surface and to repeat signal quality from day 
to day. With the fairly short-focal-length optical system and moderately large beam 
expansion, significant "speckle" was present in the collected light. The speckle 
pattern, the well-known result of diffuse reflection near the focal plane (ref. 20), 
would often overwhelm the fringe signal with its intensity. 

In view of these difficulties, the metal-coated glass inserts shown in figure 7 
were adopted. The test surface was coated with pure chromium by using vapor deposi- 
tion; a reflectance of 0.2 to 0.4 was achieved, measured with a He-Ne laser beam at 
normal incidence in air. This gave a fringe visibility in the range 
r = 0.5 - 0.7. Chromium, which has also been used by Tanner (ref. 12), was selected 
because of its excellent abrasion resistance — one coated insert has been cleaned 
over 200 times with no detectable degradation in surface reflectance. The back 
surface was treated with an antireflective coating of MgF^, which was cured at an 
elevated temperature. Figure 8 contrasts the best signals that were obtained from 
the polished aluminum inserts to those that have typically been obtained with the 
metal-coated glass inserts. 

Pure sapphire would also be expected to give reasonably good fringe visibility 
because of its high refractive index. Also, several commercial optical-coating 
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shops offer highly abrasion-resistant dielectric coatings which can be tailored to 
give the correct reflectance in oil, although the cost of this alternative is about 
five times that of chromium coating. For large surfaces, vapor-deposition is not 
practical because of the limited size of the coating apparatus. However, it may be 
possible to use a polymer film if its surface is treated appropriately. Uncoated 
optical glasses generally have refractive indices very near those of silicone oils 
(silicone oils have 1.4), so fringe visibility on these surfaces is quite low. If 
polished metal surfaces must be used, collection of the diffuse reflection must be 
minimized, e.g., by the use of a small collection aperture with less beam expansion 
and a longer focal length. Even with these precautions, erratic fringe visibility 
and subsequent higher data-rejection rates must be anticipated when metal test 
surfaces are employed. 


Vibration of the Instrument and Test Surface 

The interferometer system is not inherently sensitive to relative vibration 
between the instrument and the test surface, because such vibrations do not disturb 
the relative phase of the two interfering waveforms. However, two important cases 
in which vibration was important were noticed — one related to the use of detector 
apertures and the other to cases with substantial diffuse reflection. 

Detector apertures were used to effectively reduce the spot size, thus giving 
improved spatial resolution. This feature proved especially helpful for cases in 
which the oil film wedge angle exceeded approximately 0.05°. This optical config- 
uration is mildly sensitive to vibration because of the relative motion of the 
expanded spot on the detector aperture, which creates a change in the total light 
intensity that is incident on the photodiode. The instrument was clamped securely 
to the tunnel support rails, which alleviated the problem. 

Whenever significant amounts of diffuse reflection were present in the col- 
lected light and the sending optics focused the incident beam on the surface, the 
interferometer was extremely sensitive to the small vibration caused by the tunnel 
motor. This effect was traced to the speckle pattern created by diffuse reflection 
from the metal surface, which changes radically with the distance from the focal 
plane to the surface plane. In fact, the interferometer essentially behaves as a 
vibration measurement apparatus (ref. 20) under these conditions. As explained 
above, this situation arises when metal test surfaces are used and is accentuated by 
short, highly expanded, focused optical systems. All measures to alleviate the 
problem — longer focal lengths, unfocusing, and smaller beam expansion — resulted in a 
sacrifice of spatial resolution. The use of the coated glass inserts avoided the 
problem of speckle noise. 


Laser Output Intensity Drift 

Most He-Ne laser manufacturers specify that the total laser-beam-output power, 
after a short warmup period, will remain constant to within ±5$. Tests of three 
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different manufacturers' units of the 1-mW range verified the accuracy of the speci- 
fication. However, the time period of the output power variation proved to be in 
the range 20 to 200 sec, which was comparable to typical periods between fringes in 
our experiments. Thus, it was not feasible to remove the small mean drift via 
bandpass filtering. Implementation of a method of tracking the mean signal level by 
finding the local maximum and minimum levels proved cumbersome to implement. 

Instead, the ''scatter” (estimated at about 5$) which this problem introduced into 
the results was figured in by taking several measurements (about five) at each 
measurement point and averaging the results. This procedure had the effect of 
reducing the uncertainty in the determination of the term [C(X 2 - C(X 1 )] from about 
5$ to about 2$ to 3$ so that the overall uncertainty in the measurements was still 
estimated to be about 5$. 


Example Results in a Turbulent Boundary Layer 

Measurements were made in the Boundary Layer Wind Tunnel operated by the Fluid 
Dynamics Research Branch at NASA Ames. The test surface in the facility is level to 
within 0.1° and the static pressure varied less than 0.5$ of the free-stream dynamic 
pressure. The test section is 80 cm wide and 300 cm long, with an adjustable height 
(20 cm nominal) for pressure-gradient control. The test-wall boundary layer was 
tripped using a 0.4-mm-diameter wire located 20 cm from the inlet. Optical access 
for the interferometer was through the transparent Lexan control walls opposite the 
test surface. After the interferometer was moved to a particular measurement loca- 
tion and securely bolted to the tunnel-support rails, it was carefully focused on 
the surface and realigned. This procedure typically took about 20 min. 

The measurement procedure was as follows : ( 1 ) the insert was removed from the 

tunnel, cleaned with both acetone and lens-cleaning solution, and then wiped and 
blown clean with canned, dry compressed air. (2) Oil was applied as was described 
in the Oil Applications section using the special applicator fixture. (3) With the 
tunnel off, the insert was installed into. the test plate and bolted securely into 
place. The insert was installed flush with the surrounding test plate to within 
0.02 mm. (4) The tunnel was turned on and brought to the test condition within 
about 15 sec. Data acquisition was initiated only after the test condition was 
established. The set-point free-stream velocity was maintained to 0.25$ and the 
temperature to within 0.2°C over the period of data acquisition. 

Data reduction was accomplished using the laboratory computer. First, the 
digitized photodetector signals were displayed on the graphics terminal, and a 
contiguous portion of the record was selected where the fringe visibility was nearly 
constant for both beams. The record selected was generally about 75$ to 90$ of the 
entire record acquired. Fringe identification and the fitting of the fringe times 
to lubrication theory was then accomplished. Plotting of the fitted data was pro- 
vided to verify that no fringes were missed, although the constraints for equa- 
tion (15) also identified such records. Records meeting the constraints for 
equations (15) above were rejected; otherwise, the results C(X) from each beam 
were combined with the other measured parameters to yield the skin-friction 
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coefficient. Measurements were repeated five times at each measurement condition 
and averaged. 

Data shown in figure 9 were obtained under the conditions summarized below. 


Oil nominal viscosity 

Data record length.. 

No. fringes per record.'.... 
Freestream dynamic pressure 

Flow temperature 

Beam spacing 

Beam position (X^ - X^).... 
Glass insert reflectance... 
Aperture diameter 


....50 cs at 25°C 

400 sec 

20 to 40 

1.8 ± 0.1 in. H 2 0 

10° to 20°C 

2.77 mm 

3 to 8 mm 

20 % to 40$ 

1.5 mm 


For comparison, skin-friction estimates were also obtained by fitting measurements 
of the boundary-layer mean-velocity profiles to the Law of the Wall-Wake 
(ref. 21). The results shown in figure 9 indicate agreement between the two tech- 
niques to within 5 %, which supports the uncertainty estimate given in the section on 
uncertainty analysis. The time required to obtain these data is comparable 
(1) because the optical system at each station must be realigned and (2) because 
five oil-film measurements must be averaged so as to mitigate the effects of laser 
power drift. However, the laser oil film interferometer is a direct measurement 
technique — not subject to any assumptions regarding local-profile similarity — and is 
therefore much more widely applicable than the Law of the Wall-Wake. 


CONCLUSIONS 


An improved skin-friction laser interferometer system has been designed and 
tested. Uncertainty estimates have been used in designing the instrument and have 
been verified. A new fitting procedure developed for data analysis yielded the 
unexpected benefit of providing an objective, easily implemented method for rejec- 
tion of "bad" data records thought to be primarily the result of dirt accumula- 
tion. The final data-reduction equation used thus far was limited to the single 
case of constant pressure with two-dimensional flow, without gradients of pressure 
or skin friction and excluding body forces (e.g., gravity). 

Further refinements to the system are being pursued. (1) An automatic oil- 
application system is being developed which will allow oil to be injected under 
conditions of flow in the tunnel. (2) A third data beam, with a third photodiode 
placed behind the beamsplitter, is being added to monitor the laser-output power 
level. This information will be used to directly measure the mean output level, 
which should eliminate the need for multiple measurements. (3) Use of a larger beam 
expander and closer beam spacing are being examined in an effort to improve spatial 
resolution. 
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APPENDIX 


Four computer routines for processing the digital photodetector signals are 
included. These are listed below along with the function of each routine. 

1. SCSMOO: Digital smoothing of the photodetector output 

2. SCZERO: Fringe identification algorithm 

3. SCCFIT: Curve fit to lubrication theory 

4. SCLOFI: Calculation of skin friction coefficient 
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SUBROUTINE SCSMOO ( IRAWDT . NRAWMX . FACTOR. MODE ) 
ACRONYM: Standard Computations - SMOOth data 


PURPOSE : 

Smooth the data array - value of MODE determines the type 
of smoothing algorithm employed. 

For 3-pt . smoothing algorithm, MODE determines the no. of 
times to loop through. 

METHOD: self explanatory - see NOTES and code 
ARGUMENTS 

ARG DIM TYPE I/O/S DESCRIPTION 

IRAWDT NRAWMX I I/O The raw data to be smoothed - the smoothed 

results are written back into the same 
raw data array. 

NRAWDT - I I Length of IRAWDT array. 

FACTOR - R I Smoothing factor: (for MODE=l) 

l.= Perfectly flat curve 
0.= No smoothing 

Values between 0 . and 1 . give varying 
degrees of smoothing. 

MODE - II Determines type of smoothing to use : 

0 = single-pole or "causal" exponential 

smoothing with parameter FACTOR. 

1 to 10 = 3-pt. smoothing formula 

COMMONS USED: None 
FILES USED: None 
ERROR HANDLING: None 
NOTES : 

(1) This routine modifies the input (raw) data array. 

(2) With MODE=0, a fixed time delay is introduced into the data, 
visisble as a rightward shift of a few sample intervals. 

(3) For MODE=l to 10, no phase distortion (time delay) accrues 
from filtering. 

(4) Successive applications of the 3-pt. algorithm simply provide 
higher-order smoothing, i.e., 3-pt. smoothing used twice is 
precisely 5-pt. smoothing. 

LOCAL VARIABLES: 

VAR DIM TYPE DESCRIPTION 

SMOOTH - R Smoothed data value - returned in IRAWDT 

EXTERNAL REFERENCES: None 

STANDARDS VIOLATIONS: None. 

ENVIRONMENT: DEC RSX-11M, FORTRAN 77 

DEVELOPMENT HISTOkY: 

DATE INITIALS DESCRIPTION 
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11 '16/83 TML Initial design 

12/01/83 RVW Added 3-pt. algorithm 

01 ■'11/84 TML Loop through 3-pt. algorithm MODE times 

AUTHORS: Ted Lichtenstein, Informatics General Corp. 

Russell V. westphal, NASA-Ames Research Center 

_********************************************************************** 


INTEGER IRAWDT ( NRAWMX ) 

IF ( MODE . GT . 0 ) GO TO 250 

C**** MODE=0 : exponential smoothing (single-pole causal filter) **** 

C * If the smoothing factor is 0. or out of range, just return: 

IF ( FACTOR. LE.O. ) GO TO 999 
IF ( FACTOR. GT.l. ) GO TO 999 

C * Do the smoothing: 

FACT1M = 1. - FACTOR 
SMOOTH = FLOAT ( IRAWDT(l) ) 

DO 200 I - 1, NRAWMX 

SMOOTH - SMOOTH * FACTOR + FLOAT ( IRAWDT(I) )*FACT1M 
IRAWDT ( I ) - INT( SMOOTH+ . 5 ) 

200 CONTINUE 
GO TO 999 

C* ’ * * MODE=l : 3-pt. smoothing algorithm **** 

250 CONTINUE 

IMAX=NRAWMX-1 

C * Loop through the 3-pt. algorithm, MODE no. of times: 

DO 400 <J~ 1 , MODE 

OLD=FLOAT( IRAWDT ( 1 ) ) 

DO 300 1=1 , IMAX 

SMOOTH = ( OLD+ FLOAT ( 2 * IRAWDT ( I ) + IRAWDT ( I + 1 ) ) ) / 4 . 0 
OLD=FLOAT( IRAWDT( I ) ) 

IRAWDT ( I ) =INT( SMOOTH+O . 5 ) 

300 CONTINUE 

SMOOTH: - ( OLD+ 3.0* FLOAT ( NRAWMX ) ) / 4 . 0 
IRAWDT ( NRAWMX )=INT( SMOOTH+O. 5) 

400 CONTINUE 

999 RETURN 
END 
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SUBROUTINE SCZERO ( IRAWDT . NRAWMX, MINFIT, MAXFIT, NINTMN, 
> AVG, RHS , RMSFCT, ZERO, NZEMAX , NZERO ) 

ACRONYM: Standard Computations - find ZEROes 


PURPOSE : 

Find mean, standard deviation of data array in the interval 
defined by MINFIT, MAXFIT. Then locate positively-sloped 
zeroes of the data array in the interval. 

METHOD : 

Average and standard deviation are found using simple 
running sums of the smoothed data array. 

Zeroes are found by bracketing regions where all the data 
lie within a specified fraction of the standard deviation 
of the mean and the slope of the curve is positive. 

Then a straight line is fit to the region so located and 
the zero-crossing is. computed from this fit. 


ARGUMENTS 

ARG 

DIM 

TYPE 

I/O/S 

DESCRIPTION 

IRAWDT 

NRAWMX 

I 

I 

The raw data 

NRAWMX 

- 

I 

I 

Length of raw data array. 

MINFIT 

- 

I 

I/O 

Minimum argument for data reduction. 

MAXFIT 


I 

I/O 

Normally an input variable , but it will be 
reset if out of range. 

Maximum 

NINTMN 

- 

I 

I 

Min. no. of points in one interval used to 

AVG 

_ 

R 

0 

determine the zero-crossing line ( >= I). 
Average of smoothed signal counts. 

RMS 

- 

R 

0 

Standard deviation of smoothed signal counts . 

RMSFCT 

- 

R 

I 

Fraction of RMS to use for testing if a given 

ZERO 

NZEMAX 

R 

0 

data point is in the region near a zero . 
Abscissa locations where smoothed input wave 

NZEMAX 

_ 

I 

I 

IRAWDT has a positive-sloped zero. 

Amount of storage for ZERO( * ) in the calling 

NZERO 

_ 

I 

0 

program, in 4 -byte array elements. 
No. of zeroes found. 


COMMONS USED: None 


FILES USED: None 
ERROR HANDLING: 

Stop if there is not enough storage for the zero-crossings. 

NOTES : 

( 1 ) The zero-finder requires reasonably smooth data and a constant 
mean-line to work properly. 

( 2 ) The routine may find one zero outside the interval on each side 
if there are data near the zero just within the interval. These 
should be valid zeroes, although the calling program may elect 
to discard them. 

EXTERNAL REFERENCES: None 

STANDARDS VIOLATIONS: None. 
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ENVIRONMENT: DEC RSZ-11M. FORTRAN 77 

DEVELOPMENT HISTORY: 

DATE INITIALS DESCRIPTION 
11/ 16/80 RVW Initial design and coding 

01/12/84 TML Modified 

10/19/84 RVW Improved checking of intervals 

AUTHOR: Russell V. Westphal, NASA-Ames Research Center 

Ted Lichtenstein, Informatics General Corp. 


REAL AVG. RMS, ZERO(NZEMAZ) 

INTEGER IRAWDT ( NRAWMZ ) 

* Make sure that MINFIT and MAZFIT are in range, by resetting any 

* out of range value : 

IF ( MINFIT. LT.l ) MINFIT = 1 

IF ( MAXFIT . GT . NRAWMZ ) MAZFIT = NRAWMZ 

* Test RMSFCT and NINTMN to be sure they're reasonable 

IF ( (RMSFCT. LT. 0.1). OR. (RMSFCT. GT. 2.)) GO TO 990 
IF ((NINTMN. LT.l). OR. (NINTMN. GT. 10)) GO TO 990 

* Average the smoothed data to find the mean line on the interval; 
also, compute the rms of the (possibly smoothed) signal 

AVG-0 - 0 
RMS =0.0 

DO 300 I -MINFIT, MAZFIT 
IIND = I - MINFIT + 1 

AVG = AVG + (FLOAT(IRAWDT(I))-AVG)/FLOAT(IIND) 

DEVIAT = FLOAT ( IRAWDT (I) ) - AVG 

RMS = RMS + ( DEVIAT ’DEVIAT-RMS ) /FLOAT( IIND) 

300 CONTINUE 

RMS = SQRT( RMS ) 

* Search the data for regions which contain positive-sloped zeroes; 
bracket these regions and locate the zeroes 

NZERO=0 
IENFL-0 
NIN =0 

C *** START OF SEARCH LOOP *** 

DO 500 I=MINFIT, MAZFIT 

C — See if point is in the region near a zero 

IF ( ABS( FLOAT ( IRAWDT ( I ) )-AVG) . LT . (RMS* RMSFCT) ) GO TO 410 
C — Now if the region has been bracketed, go to fitting algorithm 
402 CONTINUE 

IF((IENFL.EQ. 1). AND. (NIN. GE. NINTMN)) GO TO 420 
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C -- Continue search for more zeroes 

405 CONTINUE 

IENFL=0 
NIN=0 
GO TO 490 

C — Loop to bracket a contiguous region near zero 
410 CONTINUE 

IF(IENFL.EQ.O) IMIN-I ! first point in interval is IMIN 
IMAX=I 

NIN=IMAX-IMIN+ 1 
IENFL=1 

C -- If this is the last point in the interval, check if 
C there's a zero in the bracketed region: 

IF ( I.EQ.MAXFIT ) GO TO 402 
GO TO 490 

C — Fit data near zero to straight line and find zero crossing 
420 CONTINUE 

C — Check that first point in this region is below the mean line , 
C -- and last point is above the mean line: 

IF ( FLOAT ( IRAWDT(IMIN) ) . GE . AVG ) GO TO 405 
IF ( FLOAT ( IRAWDT(IMAX) ) . LE . AVG ) GO TO 405 

SX=0 , 0 
SY=0 . 0 
SXX=0.0 
SXY-0.0 

DO 450 <JJ=IMIN , IMAX 

RJJ = FLOAT (JJ- IMIN) 

SX=SX+RJJ 
SXX=SXX+RJJ*RJJ 
S Y = SY + FLOAT ( IRAWDT( JJ ) ) -AVG 
SXY= SXY+ ( FLOAT C IRAWDT ( J J ) ) -AVG ) * RJJ 
450 CONTINUE 

D=FLOAT(NIN) * SXX-SX* SX 
A=(SXX*SY-SX*SXY)/D 
B=( FLOAT ( NIN ) * SXY- SX * SY ) / D 
ZDUM=-A/' B+ FLOAT ( IMIN ) 

C -- Check to be sure the zero is positively sloped, and is 
C in the interval [IMIN, IMAX] 

WRITE(4, 3000) B , IMIN , IMAX , ZDUM 
3000 FORMAT ( ' B = '.E15.7,' IMIN = ',15,' IMAX = ',15, 

> ' ZERO : ' , E15 . 7) 

IFCB.LE.O. ) GO TO 405 

IF( ( ZDUM . LT . FLOAT ( IMIN) ) . OR . ( ZDUM . GT . FLOAT ( IMAX) ) ) 

> GO TO 405 

C — If so, add it to the list of positive-sloped zeroes found 
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C -- Stop if we have run out of storage space: 

IF ( NZERO . GT . NZEMAX ) STOP ' SCZERO : LACK OF STORAGE' 

ZERO ( NZERO ) - - A/ B+FLOAT ( IMIN ) 

GO TO 405 

490 CONTINUE 
500 CONTINUE 

C *** END OF SEARCH LOOP *** 

GO TO 999 

990 STOP 'SCZERO : BAD ARGUMENT VALUE' 

999 RETURN 
END 
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SUBROUTINE SCCFIT C ZEROP, NZEROP . MODE, CFNFP . ASLFP , 

TSTFP , ERROR ) 

ACRONYM: Standard Computations - Curve FIT zero locations 


PURPOSE : 

Fit locations of a particular phase of a periodic waveform 
to the LOFI data reduciton equation. 

METHOD : 

Fit of the LOFI data reduction equation; 

ASLFP ^average of (CFNFP-I)*(TSTFP+ZEROP(I)) for 1=1, NZEROP 

The problem is to minimize the rms error of the computed 
fringe-time product with variable initial fringe number 
and base time. A grid least-square routine (GRIDLS) is 
used to perform the minimization. 

ARGUMENTS : 

ARG DIM 

ZEROP 200 


NZEROP 

MODE 


CFNFP 
ASLFP 
TSTFP 
ERROR 

COMMONS USED: 

/FITCOM/ Used looally to pass parameters to optimization routine 

FILES USED: None 


TYPE 1/0/ S DESCRIPTION 

R I Abscissa locations where smoothed input wave 
IRAWDT has a positive-sloped zero, ordered 
in time (high to low fringe numbers) 

I I No. of zeroes in data array 
I I Used to specify whether initial time base 
is to be fit or input : 

M0DE=1 : TSTFP is input 
MODE =2 : TSTFP found from fit procedure 
R 0 Fringe number offset , 

R 0 Product of time and fringe number (ANSWER) 

R 0 Time base (from fit or input) 

R 0 RMS error of fit 


ERROR HANDLING: None 


NOTES: 

(1) A simple grid-search minimization routine is used for the 
2-parameter minimization problem. 


LOCAL VARIABLES: 



VAR 

DIM 

TYPE 

DESCRIPTION 

ABSERR 

- 

R 

Absolute error convergence criteria. 

TOL 

— 

R 

Convergence oriteria on error variation 
between successive calls to GRIDLS 

ITMAX 

- 

I 

Maximum no. of calls to GRIDLS allowed 

ERRMAX 

- 

R 

Allowable error if ITMAX is equalled 


EXTERNAL REFERENCES: 

NAME DESCRIPTION AND SOURCE 

FITERR Function called by GRIDLS to evaluate error of fit 

GRIDLS Grid-search function minimization routine 
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STANDARDS VIOLATIONS: None. 

ENVIRONMENT: DEC PDP-11/RSX-11M and VAX/VMS, FORTRAN 77. 

DEVELOPMENT HISTORY: 

DATE INITIALS DESCRIPTION 
11/30/83 RVW Initial design 

3/20/84 RVW Implemented improved algorithm and numerics 

10/19/84 RVW Augmented convergence criteria 

AUTHOR: Russell V. westphal. NASA-Ames Research Center 


REAL ZEROP(200) 

REAL P ARAM ( 2 ) , DELTA ( 2 ) 

DATA ABSERR/0 . 0025/ , TOL/O . 00005/ , ERRMAX/0.05/ 

DATA ITMAX/100/ 

COMMON/FITCOM/ CFNFIT , ASLFIT , TSTFIT , ZERO( 200) , NZERO , RMSERR 
EXTERNAL FITERR 


NZERO=NZEROP 
DO 5 1=1, NZERO 
ZERO( I ) =ZEROP( I ) 


* Initial guesses for grid optimizer along with step sizes for search 

* NOTE: its easy to guess the initial fringe number, but quite difficult 

to estimate the basis time (MODE=2). The values selected will 
only work well if the basis time is near zero, i.e., if the 
fringe record begins approximately at the start of the oil flow. 


CFNFIT=FLOAT ( NZERO ) 

IF (MODE . E<? . 1 ) TSTFIT=TSTFP 
IF (MODE . EQ . 2) TSTFIT=0.0 
DELTA( 1 )=1 . 0 
DELTA (2 ) =ZERO( 2 ) -ZERO( 1 ) 


! lower bound on initial fringe number 
! time base is input (MODE=l) 

! guess for time base (MODE=2) 

! initial search grid size is one fringe 
! set time grid size to least samples 
! between zeroes 


FARAM( 1 ) =CFNFIT 
PARAM( 2 ) =TSTFIT 
OLDERR =FITERR ( MODE , PARAM) 
ITNO = 0 


5 ITNO = ITNO + 1 

CALL GRIDLS (MODE, PARAM, DELTA, DUM, FITERR) 

WRITE (5, 2000) ITNO, DUM 

000 FORMAT ( ' Iteration no.= 15, ' Error = ', F12.5 ) 

CFNFIT= PARAM ( 1 ) 

TSTFIT=PARAM( 2 ) 

* Convergence criteria: 

++ Error is less than a value ABSERR (acceptable error criteria) ; 

OR 

++ Successive values of the error are the same within a small 
fraction TOL of each other; 

OR 

++ The routine has not converged after a large number of calls 

to the grid-search routine, but the error is not terribly large. 

* IF <no convergence as defined* THEN STOP 
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NOTE: GRIDLS automatically adjusts grid search size for 
successive calls. 

IF ( DUM . LT . ABSERR ) GO TO 20 
IF( ABS(DUM-OLDERR) . LT . TOL) GO TO 20 
IF ( ( ITNO . GT . ITMAX ) . AND . ( DUM . LT . ERRMAX ) ) GO TO 20 
IF(ITNO.GT.ITMAX) GO TO 990 
OLDERR=DUM 
GO TO 15 

C * Set values for return to calling program 

20 CFNFP=CFNFIT 
ASLFP=ASLFIT 
TSTFP=TSTFIT 
ERROR =DUM 
GO TO 999 

C * Error return 

990 CONTINUE 

STOP ' SCCFIT : NO CONVERGENCE' 

999 RETURN 
END 

REAL FUNCTION F I TERR ( MODE , PAR AM) 

REAL PARAM(2) 

COMMON / FITCOM/ CFNFIT , ASLFIT . TSTFIT , ZERO( 200 ) . NZERO , RMSERR 

C * Used to compute relative rms error of the fringe-time product 
C from the average value of this quantity using input values 
C of initial fringe number PARAM(l) and time base PARAM(2). 

CFNFIT=PARAM( 1 ) 

IFCMODE.EQ.2) TSTFIT=PARAM(2) 

C * Compute rms error of the fit 

SC=0. 

SCC=0 . 

DO 1 1=1, NZERO 

C Fringe-time product for ith zero 

Cl = ( CFNFIT-FLOAT ( I) ) * ( ZERO( I ) +TSTFIT ) 

C Running-sum averages for mean and std. deviation 

SC = SC + ( Cl - SC ) / FLOAT ( I) 

SCC = SCC + (CI*CI- SCC ) / FLOAT(I) 

1 CONTINUE 

C * RMS error of fit is std. deviation of fringe-time product 
C coputed for each zero crossing 

RMSERR=SQRT( SCC-SC* SC) 

C * Average fringe-time product is taken as the answer 
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ASLFIT-SC 


C * The RELATIVE error is the quantity which is being minimized 

FITERR-RMSERR/ASLFIT 

RETURN 

END 


SUBROUTINE GRIDLS 

PURPOSE : Grid search for multi-dimensional minimization 
problems such as least-square fitting. 

SOURCE : P. R. Bevington, "Data Reduction and Analysis for 

the Physical Sciences", McGraw-Hill, pp. 212-213; 

modified by R. V. Westphal, 12/83 

USAGE : 

CALL GRIDLS ( NTERMS , A, DELTAA , CHISQR, FCHISQ ) 

NTERMS No. of parameters 

A Parameter array 

DELTAA, Increments for parameters 

CHISQR Error of fit 

FCHISQ Function to minimize 

FUNCTIONS CALLED: 

FUNCTION FCHISQ ( NTERMS, A ) -- evaluates fit error 

NOTES : 

The calling program deals with convergence of the search, 
typically by repeatedly calling GRIDLS and testing whether 
CHISQR has attained a constant value (within some tolerance) 
for two successive iterations. Note that the ever-smaller 
values for the grid steps are automatically provided by the 
subroutine on return to the caller (see statement 84). 


SUBROUTINE GRIDLS ( NTERMS , A, DELTAA , CHISQR , FCHISQ) 
DIMENSION A (NTERMS) .DELTAA (NTERMS) 

CHISQR-O. 

DO 90 J-l, NTERMS 

CHI SQ1 -FCHISQ (NTERMS , A) 

FN=0 . 

DELTA-DELTAA ( J ) 

41 A(J)=A(J)+DELTA 

CHI SQ2 -FCHI SQ ( NTERM S , A ) 

IF ( CHI SQ1 -CHISQ2 ) 51,41,61 
51 DELTA— DELTA 

A(J)-A(J)+DELTA 
SAVE=CHISQ1 
CHISQ1-CHISQ2 
CHI SQ2= SAVE 
61 FN-PN+1. 

A(J)=A(J)+DELTA - 
CHI SQ3 -FCHISQ (NTERMS , A) 
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71 

81 

84 

90 

C * 


IF(CHISQ3-CHISQ2) 71.81.81 
CHISQ1 =CHISQ2 
CHISQ2=CHISQ3 
GO TO 61 

DELTA=DELTA * ( 1 . / ( 1 . + ( CHI SQ1 -CHISQ2 ) / ( CHI SQ3 -CHI SQ2 ) ) + 0 . 5 ) 

A(J)=A(J) -DELTA 

DELTAA ( J ) =DELTAA ( J ) * FN / 3 . 

CONTINUE 

Minimum has been located - evaluate function and return 

CHI SQR=FCHISQ ( NTERMS , A ) 

RETURN 

END 
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SUBROUTINE SCLOFI ( TCOIL, TCNOM, MUNOM, NOIL, WLLSER, 

> ANGIN, DXBEAM, DELNT, TAUBAR. IER ) 

ACRONYM: Standard Computations - Laser Oil Film Interferometer 


PURPOSE : 

Implement lubrication theory data reduction equation for 
laser oil film interferometry measurement of skin friction. 
For use with Dow-Corning 200 silicone oil and applicable 
to both single- and two-beam interferometers. 

METHOD : 

Simply compute the skin friction from the linear equation. 
Oil viscosity is computed from the Dow-Corning exponential 
relation, and the effective incidence angle is computed 
from the equation for geometric optics. 


ARGUMENTS : 
ARG 

DIM 

TYPE 

I/O/S 

DESCRIPTION 

TCOIL 

- 

R 

I 

Oil temperature, in C. 

TCNOM 

- 

R 

I 

Nominal oil temperature for D-C 

MUNOM 

_ 

R 

I 

equation, in C. 

Nominal oil viscosity for D-C 

NOIL 

_ 

R 

I 

equation, in Kg/(m~s). 
Oil index of refraction 

WLLSER 


R 

I 

NOIL = 1.4 nominally for D-C oil. 
Laser wavelength in m. 

ANGIN 

_ 

R 

I 

WLLSER = 0.6328 X 10~-6 for He-Ne. 
Laser incidence angle . in deg . 

DXBEAM 

_ 

R 

I 

ANGIN = 0. for normal incidence. 

Beam separation in m. For single-beam 

DELNT 


R 

I 

LOFI , this is the distance from the 
oil film leading edge. 

1/ (difference in fringe-time product) 

TAUBAR 


R 

0 

in units of l/(fringe-s) . 
The shear stress, in N/m~2 

IER 

- 

I 

0 

Error return code. 


COMMONS USED: None 
FILES USED: None 
ERROR HANDLING: 

Check range of inputs. IER=1 means no obvious input error. 

NOTES : 

(1) This routine will only work for Dow-Corning oils! 


LOCAL VARIABLES: 
VAR DIM 

MUOIL 
ANGOIL 
OPLOIL 


TYPE DESCRIPTION 

R Oil viscosity from D-C equation, in N-s/m"2 
R Effective beam angle through oil, in rad. 

R Optical path length in the oil, in m/ fringe 


EXTERNAL REFERENCES: None 


STANDARDS VIOLATIONS: None. 
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ENVIRONMENT: DEC VMS. FORTRAN 7? 

DEVELOPMENT HISTORY: 

DATE INITIALS DESCRIPTION 
1/16/85 RVW Initial design and coding 

AUTHOR: Russell V. Westphal, NASA-Ames Research Center 


REAL PI 

PARAMETER ( PI -3 . 1415926) 

INTEGER IER 

REAL TCOIL , TCNOM , MUNOM, NOIL, ANGIN, DXBEAM, 

> DELNT , TAUBAR 

REAL MUOIL, ANGOIL , OPLOIL 

C * Initialize 

IER - 1 

C * Range check inputs: 

IF( (TCOIL. LT.O. ). OR. (TCOIL. GT. 100.)) GO TO 990 
IF( (TCNOM.LT. -200. ). OR. (TCNOM. GT. 200.)) GO TO 990 
IF( (MUNOM.LT. 0. ). OR. (MUNOM. GT.1.E6)) GO TO 990 
IF( (WLLSER . LT . 0 . ) . OR . ( WLLSER . GT . 10 . E-6) ) GO TO 990 
IF ( ( NOIL . LT . 0 . ) . OR . ( NOIL . GT . 5 . ) ) GO TO 990 
IF( ( ANGIN. LT. 0. ). OR. ( ANGIN. GT. 90.)) GO TO 990 
IF((DXBEAM.LT.O.). OR. (DXBEAM. GT.l.)) GO TO 990 
IF( (DELNT.LT. -1. E-6). OR. (DELNT. GT.1.E6)) GO TO 990 

C * Viscosity from Dow-Corning formula — kg/(m-s) - N-s/nT2 

MUOIL = MUNOM / (EXP(0.0146 * (TCOIL - TCNOM))) 

C * Incidence angle — rad 

ANGOIL = ASIN( SIN( PI ‘ANGIN/ 180. ) /NOIL) 

C * Optical path length in the oil — m per fringe 

OPLOIL - WLLSER / (2. * NOIL * COS (ANGOIL)) 

C * Linear lubrication theory - equation for shear stress — N/nT2 

TAUBAR « MUOIL * DXBEAM * DELNT / OPLOIL 

GO TO 999 

990 IER = 0 

999 RETURN 

END 
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TABLE 1.- UNCERTAINTY ANALYSIS, TERM BY TERM BREAKDOWN 

V X 2 “ X 1 } 

Data reduction equation: t = ~ r ' c(x ) - ~ c(X ' )] 

o L 2 r 


4 2 

-L «T. 


Linear, independent undertainty model: — = —• 


est sT^/T^ 


Conditions 


1 J * 0 


Flow temperature measured, manufac- 
turer equation verified 


<0.005 


Near-normal incidence measured n,. 


3 X 2 - Xl 


<0.02* c 


Beam spacing measured directly at 
focal plane 


4 C(X 2 ) - C(X.,) f 


' ! iS 

(1 -C x ) : 


C known to high accuracy 


^Spacing can be measured to ±0.05 mm -*• 2% % X 2 - X^ = 2.77 mm, the 
smallest spacing used. 

b 5C/C is the uncertainty in the determination of fringe-time product 
for a particular beam. 5C/C - 0.005 under idealized conditions. 
c C x = X 2 / Xl = C( X 2 ) /C( X i ) from equation (5). Note that 

x i - x o » i-r^ <x 2 - X 1 ) 
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(a) 



Figure 1.- Two dimensional oil film, (a) Coordinates and parameters for the film, 
(b) Laser interferometric measurement of oil film thickness. 


PHOTOCELL 



Figure 2.- Tanner’s system (fig. 2 from ref. 12). 


3 ^ 



35 


(?) POLARIZATION BEAM SPLITTER 

( 5 ) FOCUSING LENS 

(10) POLARIZER 

(R) 6328 A FILTER 

(f§) SILICON PHOTODIODE 

© MODEL SURFACE 




(a) Monson's optical system (fig. 1 
(b) Typical signal from Monson's 
ata-reduction equations. 



TO COMPUTER A/D 



APPROX. PART 

SCALE , NO. 

10 cm 1 

2 

3 

4 

5 

6 

7 

8 
M 


DESCRIPTION QTY 

0.8 mW He-Ne LASER 1 

20X OBJECTIVE 1 

f150 mm ACHROMAT LENS 1 

BEAMSPLITTER FIXTURE 1 

COLLIMATING LENS f450 1 

FOCUSING LENS f250 1 

DETECTOR ASS'Y 2 

PHOTODIODE AMPLIFIERS 2 

MIRRORED SURFACES 8 



Figure 4.- Current optical hardware. 
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GRAPHICS PRINTER 



2 EA. 

10 MB 

HARD DISK 


OUTPUT FROM 

PHOTODIODE 

AMPLIFIERS 


Figure 5.- Data-acquisition system. 
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P, arbitrary units 


0 


200 

TIME, sec 


400 



(a) Before processing. Note that the signal quality on each beam is equal, as 

provided by this design. 

Figure 6.- Typical digitized photodetector signal. 
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, sec 


ing fringe locations as identified by 
! SCZERO. 

Continued. 





(c) Fit of fringe-time history to lubrication theory, from fringes identified in 

figure 6b, using subroutine SCCFIT. 

Figure 6.- Concluded. 
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Figure 7.- Speci al oil applicator fixture and metal-coated glass-surface insert. 
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P, arbitrary units 





eOT x 'o 


□ >011 Film Interferometer 
O* Lew— of-the- Wall fit 



X - cm 

Figure 9.- Skin friction in constant-pressure turbulent boundary layer. 
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